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1  Summary 

With  this  support  we  developed,  implemented,  and  validated  a  suite  of  computational  tools  for  the 
optimal  design  of  phononic  materials  using  topology  optimization.  Specifically  we  have: 

•  Developed  an  analysis  and  optimization  code  to  handle  large  sized  3D  elastic  wave 
propagation  problems  in  heterogeneous  media  by  implementing  a  parallel  solver  to  handle 
large  systems  described  by  complex-valued  systems  of  equations. 

•  Verified  our  ideas  and  tested  the  code  on  a  number  of  numerical  examples  in  2D  and  3D 

•  Studied  many  problems  regarding  2D  and  3D  bulk  and  surface  waves  and  their  use  to 
create  wave  guiding  devices. 

•  Developed  and  implemented  an  approach  to  design  phononic  bandgap  materials  on  a 
periodic  cell. 

•  Used  our  computational  tools  to  study  the  potential  of  elastic  wave  guides  as  logic 
devices. 


The  results  have  been/will  be  presented  at. 

•  2005  ASME  International  Mechanical  Engineering  Congress  and  Exposition  (held  on 
November  2005,  in  Orlando,  FL). 

•  2006  International  Conference  on  Composites/Nano  Engineering  (Jul.  2006,  Broomfield, 
CO). 

•  2006  AIAA/ISSMO  Multidisciplinary  Analysis  and  Optimization  Conference  (Sept.  2006, 
Portsmouth,  VA) 

•  2006  Graduate  Engineering  Annual  Research  Symposium  (Mar.  2006,  Boulder,  CO). 
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•  2006  ASME  International  Mechanical  Engineering  Congress  and  Exposition  (November, 
2006,  in  Chicago,  IL). 

•  7th  World  Congress  on  Structural  and  Multidisciplinary  Optimization,  May  21-25,  2007, 
Seoul,  Korea. 


2  Objectives 

While  the  quasi-static  behavior  of  materials  and  devices  is  of  vital  importance,  the  dynamic 
behavior  can  be  even  more  significant  for  microsystems  applications  For  example,  design  goals 
for  RF  systems  are  focused  on  tailoring  the  resonance  spectrum.  Inexpensive,  batch-fabricated 
chemical  and  biological  sensors  can  be  enabled  by  micromachined  wavelength-selective 
emitter/bolometer  arrays,  but  these  must  be  precisely  thermally  isolated  (heat  is  conducted  by 
phonons).  Micromachined  acoustic  sensors,  for  example,  for  damage  detection  in  a  high- 
temperature  turbine,  require  high-resolution  incoming  waves.  Vibration  isolation  and  shock 
protection  of  microsystems  pose  significant  and  largely  unsolved  problems.  Tailoring  of  dispersion 
properties  of  materials,  in  particular  the  design  of  phononic  band  gap  materials,  and  the  design  of 
structures  with  specific  wave  guiding  and  absorption  features  has  recently  attracted  substantial 
interest.  Concerning  the  design  of  materials,  the  bulk  of  the  work  has  focused  so  far  on  theoretical 
and  experimental  studies  on  the  physics  of  band  gaps  in  periodic  materials  (Sigalas  and 
Economou,  1992;  Vasseur  et  al.  1998;  Liu  et  al.,  2000).  It  was  shown  that  the  transparency  of  a 
material  for  elastic  waves  can  be  altered  by  augmenting  a  soft  matrix  material  by  stiffer  inclusions 
Only  little  work,  however,  has  been  done  on  systematic  design  approaches  to  design  band  gap 
materials  Cox  and  Dobsen  (1999,  2000)  first  exploited  the  potential  of  topology  optimization  for 
the  design  of  photonic  band-gap  materials.  For  elastic  waves,  Sigmund  and  Jensen  (2003)  and 
Hussein  et  al.,  (2003;  2004)  have  studied  the  optimum  topology  of  unit  cells  in  infinite  periodic 
media.  The  same  research  groups  have  also  studied  waveguides  and  filters  for  bulk  in-plane 
elastic  and  acoustic  waves,  using  periodic  microstructures.  Hussein  et  al.,  (2003;  2004)  created 
defects  in  penodic  materials  to  guide  waves,  while  Sigmund  and  Jensen  (2003)  developed  a 
topology  optimization  approach  to  accomplish  this. 

Our  objective  is  to  develop  an  approach  to  design  the  layout  of  multiple  constituent  materials  that 
are  distributed  in  distinct,  often  non-intuitive,  complex  patterns  to  create  materials  and  devices  with 
tailored  phononic  properties  We  generate  these  patterns  using  topology  optimization,  building  on 
our  previous  and  ongoing  efforts  to  develop  a  general  purpose  analysis  and  design  tools. 


3  Status  of  Effort 

We  made  excellent  progress  during  this  effort.  Our  main  efforts  are  along  two  lines:  i) 
computational  method  development  and  implementation,  and  ii)  applications.  During  this  research 
we* 

•  Developed,  implemented,  and  validated  a  topology  optimization  approach  for  2D  and  3D 
wave  propagation  problems,  including  a  version  with  a  parallel  solver. 

•  Developed,  implemented,  and  validated  a  topology  optimization  approach  for  2D  phononic 
bandgap  matenals 

•  Applied  the  methodology  to  numerous  design  problems  involving  2D  and  3D  bulk  and 
surface  waves,  and  2D  phononic  bandgap  materials. 
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•  Applied  the  methodology  to  the  design  of  phononic  wave  devices  to  function  as  logic 
devices 

•  Developed  an  analysis  and  design  approach  to  study  the  effect  of  active  deformation  on 
wave  propagation  characteristics  in  phononic  bandgap  matenals  and  wave  guiding 
devices. 


4  Accomplishments/New  Findings 


Toooloav  Optimization  Approach 

In  topology  optimization  of  phononic  materials,  the  computational  domain  is  discretized  into  a  grid 
of  finite  elements  on  which  the  elastic  wave  equation  is  solved  The  general  optimization  problem 
which  is  specialized  for  the  harmonic  wave  equation  is  as  follows; 

min  r(5-,u(5.))  0<s{<\ 

si 

sJ.  ([K]+/®[c]-ft>2  [M]Ju-f  =  0  ^ 

gy-(-s/(u(j;.))<0  j  =  1,2,3... 

hk(sj,\i(sj))  =  0  k  =  1,2,3... 


In  eqs.  (1 ),  z  is  the  objective  function,  u  is  vector  of  displacements  that  is  the  solution  of  the  finite 
element  problem,  $  are  design  variables  which  define  the  material  properties  for  each  element,  [K], 
[C],  and  [M]  are  the  stiffness,  damping,  and  mass  matrices,  respectively,  co  is  the  driving  frequency 
for  the  forcing  vector  f,  and  gj  and  hk  are  constraints  on  the  problem.  The  objective  function  can 
take  many  different  forms,  but  for  the  phononic  problems,  the  following  formulation  with  its 
analytically  derived  sensitivities  seems  to  work  well: 


z  =  (j uy  L  u 

2Re([z.]wf 


-1 


-M 


3S; 


where* 


(2) 


[k]=[a:]+/®[c]-«2[a/] 

~  (3) 

a[x]  gMg_+Qj£k  m2^j\e 

ds .  ds .  ds .  ds . 

ill  i 

This  objective  function  maximizes  the  norm  of  the  displacements  at  degrees  of  freedom 
determined  by  L  This  formulation  results  an  optimization  problem  for  which  the  sensitivities  are 
easy  to  compute. 

In  order  to  determine  the  layout  of  different  materials  in  the  domain,  the  material  properties  for  the  i- 
the  element  (elastic  constants,  density)  are  varied  continuously  as  a  function  of  the  design 
variables  S\.  We  use  a  simple  linear  vanation  between  the  two  materials  of  interest,  but  it  is 
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straightforward  to  use  a  penalty  formulation  or  other  similar  techniques.  In  each  iteration  of  the 
optimization  algorithm  a  new  value  for  each  s,  is  determined  based  on  its  value  and  the 
corresponding  sensitivity  At  the  end  of  each  iteration  the  objective  function  will  be  improved  and 
eventually  the  design  variables  will  reach  or  approach  the  box  constraints.  When  all  or  most  of  the 
variables  reach  the  box  constraints  the  optimum  phononic  material  is  achieved.  We  solve  the 
optimization  problems  using  sequential  convex  programming  (SCP)  with  the  method  of  moving 
asymptotes  (MMA).  No  constraints,  other  than  the  box  constraints,  exist  on  the  problems 


Para/te/  Implementation  -  Domain  Decomposition  Approach 

In  order  to  efficiently  perform  repeated  3D  finite  element  analyses  of  surface  wave  problems  arising 
at  every  optimization  iteration  we  employ  a  parallel  solver.  The  basic  idea  behind  solving  linear 
systems  arising  from  finite  element  discretizations  of  partial  differential  equations  on  modem 
distnbuted  memory  parallel  computers  is  “divide  and  conquer.”  The  elastic  solid,  occupying  the 
computational  domain,  is  divided  into  non-overlapping  regions  (see  Figure  1): 

n  =  {Jnk 

k 


Fig.  1  Domain  decomposition  in  two  and  three  dimensions. 


Such  a  decomposition  naturally  induces  a  subdivision  of  the  finite  element  mesh.  Thus,  all  finite 
element  quantities,  such  as  the  solution  and  the  load  vectors,  as  well  as  the  stiffness,  mass,  and 
damping  matnces  may  be  restricted  to  each  subregion,  or  subdomain: 

KtJ  =  \CJ(prA(x)V(p)  )dx=  j ,  A(x)V q>j )dx ,  where  fik  2supp(p,)(Jsupp(py) 
n  n„ 

Therefore,  using  the  domain  decomposition  technique  we  transform  the  wave  propagation  problem 
in  the  whole  solid  into  wave  propagation  sub-problems  in  each  subdomain,  coupled  only  through 
the  constraints  imposing  continuity  of  the  global  solution  across  the  subdomain  interface. 

[B]u  =  0 
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These  constraints  are  imposed  weakly,  through  the  Lagrange  multiplier  technique,  so  we  have  to 
solve  the  following  linear  system: 

([ K ]  +  ico[C]  -  co2[M])u  +  [Bf  A  =  f 
[B]u  =  0 

Now  the  primal  variables  u  may  be  completely  excluded  from  the  system  above  by  solving  the 
wave  propagation  problem  locally  on  each  subdomain  in  parallel,  and  then  the  resulting  system  of 

dual  equations  in  X  is  solved  using  an  efficiently  preconditioned  parallel  iterative  solver. 

At  the  level  of  the  optimization  algorithm,  we  chose  not  to  parallelize  the  code  owing  to  the  special 
structure  of  the  topology  optimization  problems,  where  the  number  of  state  variables 
(displacements)  grows  much  faster  (as  nA3  for  a  cubical  domain  of  size  n)  than  the  number  of 
design  variables  (which  is  nA2  for  a  cubical  domain  of  size  n).  Thus  we  take  a  “master  node” 
approach,  where  only  one  selected  node  in  a  computational  cluster  is  concerned  with  optimization, 
and  contains  all  design  variables.  Thus,  the  algorithm  works  as  follows: 

1 .  Choose  an  initial  design  Xo,  set  k=0. 

2.  Send  the  material  properties  based  on  x*  from  the  master  node  to  the  rest  of  the 
computational  nodes;  compute  the  displacements  u(xi<)  a  domain  decomposition  algorithm 
outlined  above  Evaluate  the  objective  function  f(x*,  u(Xk))  and  its  gradient  cff(Xk,  u(xk))  using 
parallel  computing  and  collect  the  values  on  the  master  node. 

3.  On  the  master  node,  compute  Xk+1  =  x*  -  a^x*,  u(Xk)), 

4.  Check  whether  convergence  is  attained;  if  not,  set  k  =  k  +1  and  go  to  step  2. 


Nonlinear  Strain  Effects  for  a  Switchable  Waveguide 

As  a  mean  of  tuning/controlling  the  wave-guide,  finite  prestraining  may  be  implemented  in 
various  ways,  including  but  not  limited  to  thermal  or  mechanical  loading,  or  electric  loading 
when  used  in  conjunction  with  piezo-electric  materials.  To  demonstrate,  we  used 
mechanically  induced  finite  pre-straining  to  illustrate  the  concept  of  designing  tunable  elastic 
wave-guides  with  predefined  functional  properties.  We  assume  that  two  scales  are  present  in 
the  non-linear  dynamics  model  describing  the  behaviour  of  the  wave-guide  in  the  deformed 
configuration.  Time-independent  loads  necessary  to  achieve  the  finite  deformation  of  the 
wave-guide  are  much  larger  than  the  harmonic  loads  used  as  an  “input  signal”  to  the  wave¬ 
guide.  So  we  assume  the  following  decomposition  of  the  traction: 

t(x,t)  =  tn'(x,t)  +  £tharm(x,t) 

where  e  >  0  is  a  small  scale  parameter.  We  also  assume  that  the  displacement  fields 
involved  in  the  problem  may  be  decomposed  in  a  similar  manner: 

u(x,t)  =  unl(x,t)  +  euharm(x,t) 

We  drop  from  the  equations  of  elasticity  all  the  terms  involving  strain  squared,  which  in  turn 
allows  us  to  bypass  solving  the  underlying  non-linear  dynamics  problem  to  find  u(x,  t),  and 

instead  we  solve  a  non-linear  static  problem  to  determine  u  (x,t)  and  then  separately  solve 
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a  linear  quasi-static  problem  to  find  the  Fourier  transformation  of  xxharm (x,f)in  the  deformed 
configuration.  Complete  details  can  be  found  in  Evgrafov  et  al,  (2008). 

Applications 

We  have  applied  the  topology  optimization  approach  to  a  number  of  fascinating  problems 
including 

•  2D  Bragg  grating 

•  2D  Waveguide  (bender) 

•  2D  Pressure  wave  to  shear  wave  mode  converter  (phononic  inverter) 

•  2D  Phononic  switch 

•  2D  Phononic  multiplexer  (sort  frequencies  or  wave  types) 

•  2D  Surface  wave  filter  with  bandgap  maximization 

•  3D  Bragg  grating 

•  3D  Surface  source  waveguide 

•  3D  Bending  waveguide 

•  2D  Bandgap  meta-materials  design 

o  bulk  meta-materials  for  device 

o  plate/beam  meta-materials  for  structural  applications 

These  problems  can  be  grouped  into  three  different  types,  namely  1)  those  that  are  a  direct 
attempt  toward  the  development  of  elastic  logic,  2)  those  affecting  the  propagation  of  surface 
waves,  and  3)  those  for  the  design  of  bandgap  meta-materials 

In  the  pursuit  of  the  design  of  elastic  logic  devices  it  has  been  necessary  to  investigate  different 
methods  to  serve  this  purpose  Two  options  include  the  on/off  method  and  a  new  idea  the  mode 
conversion  method.  The  on/off  method  follows  a  similar  route  as  that  for  electronic  circuits  in  that 
the  state  depends  on  whether  the  signal  is  present  or  not.  This  resulted  in  the  design  of  a  2D 
phononic  switch  in  which  depending  on  the  state  of  two  input  ports  the  elastic  wave  would  exit  one 
of  two  output  ports.  The  second  method  is  based  on  an  idea  that  the  mode  of  the  wave  (pressure 
or  shear  wave)  acts  as  the  on  or  off  state.  This  idea  resulted  in  the  design  of  a  2D  phononic 
inverter  which  converts  one  wave  type  into  another  and  the  design  of  a  phononic  multiplexer  which 
sorts  incoming  waves  by  type  into  different  output  ports.  Both  methods  have  shown  promise  in 
reaching  the  larger  goal  of  the  design  of  elastic  logic. 

Surface  wave  device  design  has  been  a  major  focus  because  of  their  applicability  to  real  device 
designs.  A  current  practical  example  is  Surface  Acoustic  Wave  filters  used  in  may  electronics 
applications  We  reproduce  a  similar  device  and  improve  on  its  basic  design  by  optimizing  the 
bandgap  at  which  it  operates.  Our  interest  extends  further  than  that,  however,  to  the  design  of  3D 
surface  waveguides  which  direct  the  propagation  of  surface  waves  in  a  similar  manner  as  the  2D 
examples  and  may  as  well  lead  to  elastic  surface  wave  logic  devices.  A  particularly  challenging 
example  has  been  the  3D  waveguide  shown  in  Fig.  2.  A  half-space  is  approximated 
computationally  with  the  use  of  wave  energy  absorbing  boundary  conditions  into  which  a  surface 
wave  propagates  in  a  predetermined  direction.  In  a  homogenous  material  the  surface  wave 
passes  through  the  computational  domain  unimpeded.  In  this  problem,  however,  it  is  desirable  to 
have  the  wave  exit  the  domain  at  a  location  perpendicular  to  the  entrance  location  thereby 
requiring  the  wave  to  bend  a  comer.  The  optimization  problem  is  set  up  such  that  the 
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displacements  are  maximized  at  an  output  location.  This  problem  has  been  solved  before  as  a 
topology  optimization  problem  by  Bendsoe  and  Sigmund  (2003)  for  bulk  waves  using  a  two- 
dimensional  model  but  has  never  been  done  for  surface  waves  in  three  dimensions.  The  results  of 
this  problem  can  be  seen  in  the  figure  which  shows  isosurfaces  located  at  the  interface  between 
the  two  design  materials  silicon  (blue)  and  aluminum  (red).  The  topology  of  the  design  is  non- 
intuitive  but  can  be  imagined  it  guides  the  surface  wave  by  reflecting  it  off  of  the  curved  surfaces  in 
some  complex  manner. 


Fig,  2  Final  design  of  a  3D  surface  wave  bending  waveguide  (left)  causes  an  incoming  surface 
wave  to  bend  a  comer  (right). 


We  also  worked  extensively  to  design  bandgap  meta-materials  through  which  elastic  waves  of 
certain  frequencies  will  not  propagate.  These  meta-materials  are  based  on  arrayed  unit  cell 
structures.  It  is  the  design  of  structure  of  the  unit  cell  that  determines  which  waves  will  propagate 
and  which  will  not.  So  far  we  have  designed  frequency  filtering  2D  bulk  meta-materials  as  well  as 
plate/beam  structures  which  serve  the  same  purpose;  an  example  is  shown  in  Fig.  3. 


A 


H  Aluminum  H  Silicon 


▲ 


Fig.  3  Design  of  a  bandgap  meta-material  designed  to  have  a  bandgap  at  a  prescribed 
frequency  range 
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Such  technology  is  currently  promising  because  of  its  application  to  materials  with  a  negative 
refraction  index  which  have  been  shown  to  allow  for  the  design  of  such  things  as  invisibility  cloaks 
or  superlenses.  The  advantage  we  have  is  that  unlike  over  current  groups  in  this  research  area  we 
have  the  ability  to  design  for  an  optimal  unit  cell. 
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7  Interactions/Transitions 

Results  of  this  research  have  been/will  be  presented  at: 

•  2005  AFOSR  Contractors  Meeting  (Santa  Fe,  NM,  Aug.  2005),  oral  presentation 

•  2005  ASME  International  Mechanical  Engineering  Congress  and  Exposition  (Orlando,  FL, 
Nov  2005). 

•  2006  ASME  International  Mechanical  Engineering  Congress  and  Exposition  (Chicago,  IL, 
Nov.  2006). 

•  Project  results  were  presented  to  and  discussed  with  the  Program  Manager,  Dr.  Clark 
Allred  at  a  visit  to  CU  in  April  2006. 

•  2006  AFOSR  Contractors  Meeting  (Seattle,  WA,  Aug.  2006),  to  be  presented  as  a  poster 

•  1 1th  AIAA/ISSMO  Multidisciplinary  Analysis  and  Optimization  Conference  (Sep  6-8,  2006) 

•  2006  International  Conference  on  Composites/Nano  Engineering  (Jul.  2006,  Broomfield, 
CO). 

•  7th  World  Congress  on  Structural  and  Multidisciplinary  Optimization,  May  21-25,  2007, 
Seoul,  Korea. 

•  ASME  IMECE  2006,  IMECE2006-15071. 

8  Discoveries,  Inventions,  or  Patent  Disclosures 

None 


9  Honors/Awards 

•  Cory  Rupp  won  best  presentation  awards  at  the  annual  Graduate  Engineering  and 
Research  Symposium  at  CU  in  2006  and  2007. 

•  Martin  L.  Dunn  was  awarded  the  Victor  Schelke  Endowed  Chair  at  CU. 

•  Kurt  Maute  was  awarded  the  College  of  Engineering  Joseph  H.  Smead  Faculty  Fellowship 
atCU. 

•  Martin  L.  Dunn  was  elected  a  Fellow  of  the  American  Society  of  Mechanical  Engineers 

•  Our  paper  “Large-Scale  Topology  Optimization  of  Phononic  Meta-Materials  for  Surface 
Wave  Devices,”  in  Proceedings  of  7th  World  Congress  on  Structural  and  Multidisciplinary 
Optimization,  May  21-25,  2007,  Seoul,  Korea  was  selected  as  Finalist  for  Best  Paper 
Award. 
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